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Abstract 



Second-order spin-wave expansions are used to compute the ground-state 



'^ ' energy and sublattice magnetizations of the quantum one-dimensional Heisen- 

Q ' berg ferrimagnet with nearest-neighbor antiferromagnetic interactions and 

two types of alternating sublattice spins Si > 82- It is found that in the 

'ijj- ' extreme quantum cases (51,5*2) = (1,1/2), (3/2,1), and (3/2,1/2), the es- 

7—i \ 

^^ ' timates for the ground-state energy and sublattice magnetizations differ less 
(N 



than 0.03% for the energy and 0.2% for the sublattice magnetizations from 
the recently published density matrix renormalization group numerical calcu- 
lations. The reported results strongly suggest that the quantum Heisenberg 
^ ' ferrimagnetic chains give another example of a low-dimensional quantum spin 

o' 

O ' system where the spin- wave approach demonstrates a surprising efficiency. 
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Recently, two scientific groups have published theoretical results concerning the model 
of the one-dimensional Heisenberg ferrimagnet (IDQHF) containing two alternating site 
spins ( Si > S2) per unitary cell and nearest-neighbor antiferromagnetic bonds [|1]-|^. The 
presented linear spin-wave analysis demonstrates a substantial reduction of the classical 
sublattice spins imposed by the quantum fluctuations. In the extreme quantum cases with 
(^i, S2) = (1, 1/2), (3/2, 1), and (3/2, 1/2), the quantum reduction of S2 was shown to be 
about 61%, 46%, and 37%, respectively. The density matrix renormalization group results 
01 point towards a smaller reduction. The cited values are respectively reduced to 42%, 36%, 
and 28%. For comparison, in the square-lattice Heisenberg antiferromagnet the discussed 
reduction is about 39%. It is clear that the linear theory overestimates the role of the zero- 
point spin fluctuations. In this respect, an open question is if the qualitative picture based 
on the LSWT (and concerning, in particular, the structure of the elementary excitations) 
can reflect the real situation at all. The purpose of the present paper is to throw some light 
on the above problem through an explicit study of the large S series for the ground state 
energy and sublattice magnetizations up to second order in 1/S. 

The Hamiltonian of the IDQHF with two spins Si > S2 per unitary cell and nearest- 
neighbor antiferromagnetic couplings reads 

H = 2^ Si„S2n+6, (1) 

n,(5 

where the intergers n number the cells and the vector 5 = ±1/2 connects the two nearest- 
neighbor spins. The size of the elementary cell and the exchgange interaction are unities. 
In what follows we frequently use the notations S1/S2 = w, S2 = S. 

We will use the Dyson-Maleev representation for the site spin operators. After some 
standard procedures 0, including the Fourier and Bogoliubov transformations and the nor- 
mal ordering of the boson operators, the spin Hamiltonian can be recasted to the following 
form 

H = Ho + V, V = ci + V2 + Vdm, (2) 



Ho = {-2wS' + 2Sco) N + 2SY^ [4"^«I«fc + 4^^PlPk] , (3) 

k 

^2 = E [Vi'^^l^k + Vt'^^% + Vt^'^aipl + vt^aA] , (4) 

k 

where Cq = -[1 - {l/N)Y.kek]{w + l)/2, Ci = -2((72 + gl) - 2gMw + l)w-'/\ Qi = 
-(l/2iV)Efc7fc%Afc, 92 = -(l/2) + (l/2iV)Efcl/efc, e^ = (1-%')'/', % = 2^,w^'^ /{w + l), 
7fc = cos{k/2), and A^ is the number of cells. 

Hq is the quadratic LSWT Hamiltonian. The boson operators ak and jSk describe two 
types of elementary excitations with energies 



4"'^) = 2^u;i"'^) = 2^ 



w + I (w — 1) 



(5) 



^2 2 

The a excitations are gapless (1^^ ~ k^ for small k) and describe magnons in the sector 
with a total spin (^i — S2)N — 1, whereas the /? excitations are gapful {luI = w — I + 0{k'^)) 
and belong to the sector (^i — 5*2 ) A^ + 1. 

The interaction V contains three different terms: The constant ci gives the first-order 
correction to the ground state energy. The quadratic interaction V2 introduces four vertex 
functions defined as follows 

Vk = —^- - ^IkTjk + -^^^k 2 -^IkVk (6) 

Cfc V V^ \/w tk\ Vw J 



T/(3) 9i (^ w + l \ g2 (w + l w-1 \ 

Vk = 27fc ^rik -^^Ik + -^^ik^k - 2r]k (7) 



tf = Vr + 29,Hi^7.. ff ' = fi" - 29,Hi^. (8) 



W \/W 



For the model under consideration the off-diagonal terms of V2 do not vanish due to the 
inequality 5*1 > 6*2. Note also that V2 is a non-Hermitian operator for the same reason. 

Vdm is the Dyson-Maleev normal ordered quartic interaction, containing nine vertex 
functions V^^^ = V^i2;34, i = 1, •••, 9. We have adopted the symmetric form used in Refs. P,|^ 
and the convention (fci, ^2, ^3, ^4) = (1, 2, 3, 4). 



Now, let us represent the series for the ground state energy and the magnetization of the 
first sublattice in the following form 



Eq 

N 



~2wS^ + 2Sco + Ci + ■% + ... 



mi = wS + bo -\ — - 



2S 



+ .... 



(9) 



(10) 



2S (25)2 

The coefficient bo = —g2 gives the spin reduction in the LSWT. The first-order correction 
for mi is related to the off-diagonal terms in V2, and a simple calculation gives the following 
result 



bi 



1 



1 



2(w + l)N 



T.ivt" + vt)% 



(11) 



The coefficients C2 and 62 will be calculated using the Rayleigh-Schro dinger perturbation 
theory. A straightforward calculation gives the following two contributions to C2 = C21 + C22 
resulting from V2 and Vdm, respectively 

1 1 T/(3)t/(4) o 1 T/(7) T/(8) 

'"""^^^iV^^^' '''--^^llPt.\i + e2 + e, + e, ^^'^ 

. Here 5^| is the Kronicker function. The symmetric vertex functions V"/2-34 and 1^12-34 read 



V, 



(7) 



12;34 



Wi234{a;i 0:4(71-3 - w ^/^a;37i) - (w^/^7i_3-4 - 2:371-4) 



+ X2 



3^4(72-3 - W ^/^X372) - (w^/%2-3-4 - 3:372-4) } 



^(8) 



^12734 = ^1234(3:2 (3:371-3 - w ^/%i) - a:4(w^/\37i_3_4 - 7i_4) 

+ Xi 



[3:372-3 - w ^/^72) - 3:4(^^/^X372-3-4 - 72-4) }, 



(13) 



(14) 



where W1234 = M1M2M3M4, Uk = {I + efe)/2efe, and xu = r]k/{l + Cfc)- 

To calculate &2 "we introduce in /J a staggered magnetic field for Si spins through the 
operator —hi Yl,n Sin- Then the required second-order correction for mi can be deduced 
from m^i = —{l/N)[dEQ {hi)/dhi]\i^ ^q, where Eq (hi) is the second-order correction to 
Eq in a finite hi. As a matter of fact, hi reduces the structure factor to rjk 1 — > % = 



27fcW^^^^/(l + 'w + hi/2S), and this is all one needs to find all necessary derivatives. In some 
cases, it is better to use an infinitesimal perturbing field hi and the related perturbed ground 
state I 0,,,) =1 0) + [hi/2S{w + 1)] Ekim/elHPl I 0) + 0(/i?). 

Using the latter approach, one can easily find two of the contributions to 62 = &21 + 
^22 + ^23 + ^24 related to the interaction Vdm 



J21 



{w 


2 

&22 


1 

ATS 
~ 7 


1-4 

2ti;V2 

^„ _L 1^ 


'^12;34 '^34;12 

(ei + 62 + es + 

1 V ^ r34 

3 An 2^ "12 .. ^ 


£4)2 


/I 

Ui 

2;34 


+ 

1 


1 
£2 

5 


+ 


1 


+ 


1 
£4 



(15) 



(16) 



W 



12:34 



v; 



(7) 



12;34 



^2 1/(5) , 73..(2) 

2 '^34;12 "T 2 '^34;12 
^2 ^3 



+ 



7l-r.(6) ,74 ..(3) 

2 ^12;34 "T 2 '^12;34 
Ci 64 



T/(8) 
•^ 34:12- 



(17) 



K7) 



The vertex functions appearing in the last equation have the same structure as V^i2;34 and 
V^2;34 and we will not present here their explicit expressions. 

For the other two contributions to 62? which are related to the interaction V2, we use the 
first of the mentioned methods. The result reads 



1 _ vt'^vt^ 



^23 



{w + lfN 



E 



^k 



'' {w + 1)2 iV V ek 



u^vlf^ 



ujf^v^y 
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where U, 

r(4) 



(3) 



-2gi'^k/tk - 25-2(1 - ll)Vk/4 - 93 [{w + l)w ^/2^fe + {w- l)w ^/^^fcCfc - 2r]k]/ek, 



r{3) 



up = Ul' + 2g^{w - l)w-^/2^fe, 53 = -{1/2N) Y.k Vk'^/(^k^- Notice that the equation mi + 
7712 = 'S'l — 5*2, which is connected to the conservation law for the total magnetization, 
is fulfilled order by order in the spin-wave series, so that it is enough to know only the 
corrections for one of the sublattice magnetizations. 

The results for the series for a number of combinations (5*1,5*2) are presented in the 
tables. Surprisingly, even in the extreme quantum cases (3/2, 1) and (1, 1/2), the deviations 
from the density matrix renormalization group results are less than 0.033% for the energy 
and 0.2% for the sublattice magnetizations. One can also notice that the insrease of the 
ratio w = 5*i/5*2 for fixed 5*2 = 5* = 1/2 leads to a rapid improvement of the series. This 



tendency is expected because for large w the quasiclassical spins on the first sublattice act 
as an effective field on the 5*2 = 1/2 spins. 

It is instructive to compare the series for IDQHF with those for the square-lattice Heisen- 
berg antiferromagnet (2DQHA) ||^. Let us take the series for w = 2 and for the second 
sublattice magnetization 1712- 

IDQHF: -^ = -452 -0.436456 X (25) -0.024384 + 0.006518 X — + ... (19) 

2DQHA : -^ = -45^ _ 0.315895 x {2S) - 0.024948 + 0.000866 x — + ... 

IDQHF : -ms = S- 0.3048865 + 0.1212303 x — - 0.0224602 x — ^ + ... 

2jS (20) 

2DQHA : m = S - 0.1966019 + x -^ + 0.00348 x — ^ + ... 

25 (2>j) 



The Eq spin- wave series for the compared models have similar structures.. The 1/S cor- 
rection in IDQHF is somewhat larger, but note that the LSWT reduction is also larger for 
the ID model. The values of the coefficients ci are very close to each other. As to the 
sublattice magnetizations, the main difference comes from the lack of 61 corrections in the 
2DQHA which is connected with the symmetry of the square lattice. Note that up to the 
first order the quantum spin reductions in the models are approximately one and the same. 
As a whole, the spin-wave series for the IDQHF demonstrate features which are compatible 
with those of the 2DQHA spin- wave series. It is well-known that the spin- wave results for 
the ground state parameters of the 5* = 1/2 square-lattice Heisenberg antiferromagnet are 
close to the most precise numerical estimates 0. The reported results strongly suggest that 
the quantum Heisenberg ferrimagnetic chains give another example of a low-dimensional 
quantum spin system where the spin-wave approach demonstrates a surprising efficiency. 
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TABLES 
TABLE L The coefficients of the spin-wave series for the ground state energy per cell 

eo = Eq/N of ferrimagnetic chains with two different spins: Si = 'WS2, S2 = S, w > 1; 

eo = -2wS^ + co{w) X (2S) +ci{w) + C2H x {2S)-^ + O [(23)'^)] 



{Si,S2) w = S1/S2 



co{w) 



ci{w) 



C2{w) 



eo (SWT) eo (DMRG) 



-0.41403688(9) -0.03950436(9) 0.00874156(4) -3.86320737 -3.86192 

-0.43645559(0) -0.02438442(5) 0.00651791(7) -1.45432210 -1.45408 

-0.45803557(4) -0.01179278(6) 0.00283359(4) -1.96699477 -1.96724 

-0.46862597(4) -0.00691029(8) 0.00139788(3) -2.47413839 

-0.49305421(4) -0.00037529(8) 0.00002099(5) -9.49340852 
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(1,^) 
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V2' 2> 
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(2,^) 
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TABLE II. The coefficients of the spin-wave series for the sublattice magnetization 
mi of ferrimagnetic chains with two different spins: Si = WS2, S2 = S, w > 1; 
mi=wS + boiw)+bi{w) x {2S)-^ + b2{w) x {2S)-^ + O [{2S)-'^)] 



{Si,S2) w = S1/S2 



bo{w) 



bi{w) 



b2{w) 



mi (SWT) mi (DMRG) @ 



-0.46005842(4) 0.22694904(9) -0.02899689(7) 1.14616688 1.14427 

-0.30488650(5) 0.12123033(0) -0.02246016(3) 0.79388366 0.79248 

-0.18644025(0) 0.05316804(7) -0.01006592(1) 1.35666188 1.35742 

-0.13512460(0) 0.03000271(9) -0.00504066(8) 1.88983745 

-0.02818572(2) 0.00152353(6) -0.00007830(4) 8.97325951 
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